SUBATECH-97-27 



Relativistic Transport Theory for Systems Containing 

Bound States 

P. Rehberg 

Institut fur Theoretische Physik, Universitdt Heidelberg, 
0^ , Philosophenweg 19, D-69120 Heidelberg, Germany 



9" 

43 



X 



and 



S ! SUB ATE CH 



Laboratoire de Physique Subatomique et des Technologies Associees 
UMR Universite de Nantes, IN2P3/CNRS, Ecole des Mines de Nantes 
4 Rue Alfred Kastler, F-44 070 Nantes Cedex 03, France * 



> 
On 
CO 

2 ' Abstract 

oo 
On 



Using a Lagrangian which contains quarks as elementary degrees of 
freedom and mesons as bound states, a transport formalism is devel- 
qj , oped, which allows for a dynamical transition from a quark plasma to 



a state, where quarks are bound into hadrons. Simultaneous transport 
equations for both particle species are derived in a systematic and con- 



H ' sistent fashion. For the mesons a formalism is used which introduces 

off-shell corrections to the off-diagonal Green functions. It is shown 
that these off-shell corrections lead to the appearance of elastic quark 
scattering processes in the collision integral. The interference of the 
processes qq — * tt and qq — > tt — * qq leads to a modification of the 
s-channel amplitude of quark-antiquark scattering. 
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I. INTRODUCTION 



The search for the quark-gluon plasma (QGP) as a new state of strongly inter- 
acting matter, where quarks and gluons become deconfined and can move freely, 
has been one of the most important subjects of nuclear and particle physics during 
the last decade. (For an overview see PJ.) While experiments to detect the QGP in 
heavy ion collisions have been performed and are still performed today, the theoret- 
ical modeling of these experiments still contains many open questions, one of these 
being the description of the transition from a gas of quarks and gluons to a gas of 
hadrons. While hydrodynamical models |2| have been successfully applied as phe- 
nomenological approaches, the derivation of equations of motion from a microscopic 
Lagrangian has not yet been put forward sufficiently in order to treat a system, 
that initially contains quarks and gluons and which subsequently transforms into a 
system of hadrons. 

Although the Schwinger-Keldysh formalism provides a convenient tool for 
the derivation of transport equations from microscopic models, there are still open 
questions to be addressed before it can be applied to a hadronizing QGP. One 
imagines that the underlying Lagrangian should, as in quantum chromodynamics 
(QCD), not contain hadrons as elementary degrees of freedom, since these should 
appear as bound states. On the other hand one has to consider both quarks and 
hadrons, since the former constitute the relevant degrees in the initial, the latter 
those in the final stage of the collision. Hadrons have thus to be described as effective 
degrees of freedom and the formalism has to allow for a dynamical transition between 
these two regimes. This has to be accomplished in a systematic and consistent way 
in order to avoid double counting. 

The second requirement for the theory is that it should be possible to derive 
collision terms for higher order processes. To be definite, say e. g. that in lowest 
order one obtains an interaction between quarks and pions via the process qq <-> n. 
In higher orders also the processes qq <-> qq || and qq <-> nn |3j give important 
contributions and should be derivable within the theory. A first sketch of this 
programme has been given in Ref. J7|. As it turns out, this is by far the harder task 
to fulfill. 

The description of nonequilibrium matter in terms of quasiparticles is in so far 
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problematic as it is strictly valid only in the low density limit. In Ref. || however, 
a formalism has been developed which to a certain extent overcomes this problem 
and has been successfully applied to semiconductor physics. The basic trick of this 
approach is a decomposition of the off-diagonal Green functions of the Schwinger- 
Keldysh formalism into a singular and a regular part. It turns out that in leading 
order of a gradient expansion the singular part obeys a Boltzmann equation, whereas 
the regular part gives corrections to the collision terms. It has been shown already 
in Ref. ||, that these corrections are essential if collision terms for higher order 
processes like those described above have to be considered. 

In the present paper, the approach of Ref. || is thus applied to the dynamical 
formation of pions out of a quark plasma. As a microscopic model for the interac- 
tion the Nambu-Jona-Lasinio (NJL) model 0|lOj is used. Although this model fails 
to provide a complete description of a QGP since it is not confining and does not 
contain gluons, it has the same internal symmetries as QCD and thus gives a good 
description of the low energy mesonic sector that is governed by strong interactions. 
The underlying principle of this approach is the concept of chiral symmetry break- 
ing, which forms also the basis of other phenomenological models such as chiral 
perturbation theory and the sigma model [JIT]], and which works as follows: In the 
limit of vanishing current quark masses, QCD is invariant under transformations of 
the form if> — > exp (— i8 a Xal 5 ) if), where 6 a is an arbitrary vector in flavor space and 
Xa the generators of the group SU(Nf). This symmetry, in turn, is dynamically 
broken and thus not observed in nature. As a consequence, due to the Goldstone 
theorem, Nj — 1 massless bosons appear, which for Nf = 2 are identified with the 
pions, and for Nf = 3 with pions, kaons and the rj. The finite mass of these particles 
observed in nature comes about due to a (small) explicit chiral symmetry breaking 
by finite current quark masses. It turns out that the concept of chiral symmetry is 
already sufficient to describe the low energy phenomena of strong interactions. 

In the NJL model ||10|| , this is implemented as follows: one starts with a La- 
grangian containing free quarks, which interact via a chirally symmetric contact 
interaction. This interaction in turn leads to a spontaneous breakdown of chiral 
symmetry. As a consequence the pions, which appear as bound states of quarks and 
antiquarks, become massless. The physical picture is thus that one has pions as 
ground state and constituent quarks with masses of about 300-350 MeV. At finite 
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temperature the situation changes, however, since chiral symmetry is restored at 
T « 200 MeV. In this case, the quark mass drops down to the current quark mass, 
while mesons become unstable resonances. Within the present accuracy of lattice 
calculations, chiral symmetry is restored at the same temperature where deconfine- 
ment happens The lattice data show as well, that the temperature behaviour 
of the meson masses is qualitatatively the same for both lattice QCD and the NJL 



model, as can be seen by comparing the data of Ref. [12] with those given in Ref. [10 



Due to this reasons, its simplicity and since the appearance of bound states 
is well studied in equilibrium, the NJL model provides an ideal starting point for 
studying the nonequilibrium behaviour of bound states. The results obtained here 
may nevertheless be of larger significance than the application to a specific model. 

The derivation of transport equations for the quark sector of the NJL model 



has already been done in Refs. |Ji3| , |14| and numerical solutions have been given in 
Refs. ||15|| . These works, however, do not include the formation of bound states. 
An approach to bound state formation has been attempted in Refs. ||L6|,r^|, but 
without addressing the problem of computing the collision terms for higher order 
processes. On the other hand, there have been several attempts at constructing 
transport theories for mesonic Lagrangians that treat mesons as point like interacting 
particles [|T^,|T^]. The present approach goes beyond these works in that it unites 
both the quark and mesonic degrees of freedom in a single model, and it is required 
to account for the fact that the mesonic states are bound states of quarks below 
the Mott temperature |2tJ and resonances above this. It is further attempted to 
compute the collision terms beyond the leading order. 

This paper is organized as follows: In Section H, the model is introduced. Since 
the NJL model is a strongly interacting theory, it has to be treated by an expansion 
in the inverse number of colors, 1/N C , rather than by an expansion in the coupling 



constant pl|j22| . Thus mesons, which appear already in the lowest order of this 
expansion, are introduced as effective degrees of freedom in the beginning. In Sec- 



tion |II A| , the equations of motion for the quark Green functions in coordinate space 
are obtained and the derivation of the Boltzmann equation for the quark densities 
in the standard formalism [fT4]] , which is sufficient for the present purpose, is briefly 
reviewed. The equations of motion for the mesonic Green functions are given in Sec- 



tion II B. These are treated further and cast into a Boltzmann form in Section III. 
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Section IV is dedicated to the explicit derivation of collision terms. In Section IV A 
the lowest order of the mesonic collision term is shown to contain processes of the 
type qq «-> tt, which are kinematically allowed at high temperatures. The lowest or- 
der collision integral for quarks is computed in Section [IV B| . Since in this collision 
term the off-shell corrections to the off-diagonal meson Green functions appear, this 
calculation is divided in two parts. First the contributions of the singular (quasi- 
particle) part of the Green functions is shown to yield processes of the type qq <-> tt. 
Afterwards the contributions of the regular part are shown to give processes of the 
type qq <-» qq and qq <-> qq. It is essential to consider the off-shell corrections in 
order to obtain these processes at all. It is, however, also shown that the form of 
the off-shell corrections given in Ref. || is insufficient in order to describe the cor- 
rect s-channel scattering amplitude of quark-antiquark scattering. Conclusions are 
presented in Section |V|. 



II. FUNDAMENTAL EQUATIONS OF MOTION 

The starting point for the following investigations is the NJL Lagrangian, which 
in its two flavor version reads HTU1 



C = vp(i$ - m )i/j + G {i>ip) 2 + (i>i"f 5 Tijj) 



(2-1) 



Here ip denotes the quark wave function, which implicitly contains color and flavor 
degrees of freedom, and f are the Pauli matrices in flavor space. The equations of 
motion for the quark fields obtained from Eq. ( |2.1| ) are 



-2G 



(ipip) + (ipi^f^i^r 4) 



ip(—i $ — m ) = -2Gvp (ipip) + (ipi>y 5 T'ip)i'y 5 T 



Sigma and pion fields can be introduced via |T0,O,[T7| 



7T 
7T 



= -2GipT a tp 

-2Gtpi-f 5 T 3 ip = -2G^r ^ 
n + ir 2 



-2Gi>i l5 



V2 
Ti - ir 2 

V2 



2G^r_^ 



(2.2a) 
(2.2b) 



(2.3a) 
(2.3b) 
(2.3c) 

(2.3d) 
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where, in order to simplify the notation in the following, also the matrices IV, T 
and r± have been defined. Note that these definitions imply that the a and tt° fields 
are real, whereas ir + and 7r~ form a complex field with 7r + * = 7r~. 

The equations of motion for the quark fields fl2.2| ) can be rewritten as 



(i $ - m )^ = [aY a + 7r°r + 7r+r_ + 7r-r + ) ^ (2.4a) 

-mo) = i) {oT a + 7T r + 7T+r_ + 7T-r + ) . (2.4b) 

In this form they can, together with the definition of the sigma and pion fields (|2.3|), 



be obtained from the Lagrange density 

C = ${i@- m )ij - ^ (a 2 + tt° 2 + 27T+7T-) (2.5) 
- $ (aY a + vr°r + 7r+r_ + 7r-r+) ^ , 



which is equivalent to the original Lagrangian ( [2.1| ). A canonical quantization of the 
Lagrangian Q2.5p fails however, since ( |2.5| ) does not contain the time derivative of the 



mesonic fields and thus inhibits the definition of canonical momenta for the meson 
fields. One possible way out of this problem is to quantize (|2.5| ) using the path 
integral approach, as is done in Ref. [|17]]. Here, however, an infinitesimal kinetic 



energy for the meson fields is added to the Lagrangian, 

1 



= £ + 



(d a) 2 + ((V ) +2(«9 7r + )(d 7r-) 



(2.6) 



quantization is done canonically and the limit e — >• is taken at the end.[] Using the 
generic notation 0& for the meson fields, one then obtains the associated canonical 
momenta 

dCl 

Pk = TTTTT—T = e<9 O 0fe , (2.7) 
O (9 O 0fc) 

which allow for the introduction of canonical commutation relations at equal time 



One could also think of adding other terms containing spatial derivatives of the fields 
to the Lagangian in order to maintain Lorentz covariance even for 6 / 0. However, this 
procedure would not change the following calculations and for the present purpose Eq. fl2.6| ) 
is sufficient. 
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M*)Mv) - Mv)Mx) = ° ( 2 - 8a ) 

Pk(x)pi(y) - pi(y)p k (x) = (2.8b) 
<p k (x)<pi(y) - pi(y)<f) k {x) = i5 k i5 3 (x - y) . (2.8c) 



Using these relations, the canonical quantization of ( |2.6| ) can proceed in a standard 
fashion ||23|| . 



A. Quark Fields 

The transport equations for quark fields in the NJL model are derived in 
Refs. [ |T3| , |T1| , |T7| . Here the main results of these works are listed for further reference 



and in order to define formalism and notation. 

The Green functions for the quark fields in the Schwinger-Keldysh formalism 
are defined by 

S+(x,y) = -i(T[iP(x)ij(y)]) (2.9a) 

S-{x,y) = -i(f[iP(x)>ij){y)]) (2.9b) 

S > (x,y) = -i(^(x)ij}(y)) (2.9c) 

S<{x,y) = i(${yW{x)) , (2.9d) 

where T and T are the time ordering and anti-time ordering operators, respectively: 

T [A{x)B{y)\ = 6(x - y )A(x)B(y) ± 9(y - x )A(y)B(x) (2.10a) 
f [A{x)B{y)) = % - x )A(x)B(y) ± 6(x - y Q )B(y)A(x) . (2.10b) 

Here the upper (lower) sign refers to bosonic (fermionic) operators. In order to unify 
Eqs. (|2.9|), the field operators are analytically continued to the complex time plane 
and the Schwinger-Keldysh contour C (see Fig. []]) is introduced ||]. This contour 
runs from the time to, where the initial conditions are specified, to an arbitrary 
time larger than any other time considered (conveniently set to +oo) and back. The 
initial time to will in the following be shifted to — oo. This can be done since one has 
two time scales involved, namely the mean collision time and the expansion time. 
An initial time to = — oo in this sense means thus a time which is large compared to 
the mean collision time, but still small compared to the time scale of the expansion. 
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For a more formal argument for this procedure the reader is referred to the paper 
of Botermans and Malfliet |4[] . 

The contour ordering operator Tc is defined as 



T c [A{x)B{y)\ = 6 c (x , y )A(x)B(y) ± 6 c (y , x )B(y)A(x) 
using the contour theta function 



(2.11) 



Oc(xo,yo) 



1 if Xq is 'later' than y with respect to C 
if x is 'sooner' than y with respect to C 



(2.12) 



Using this definition, Eqs. 



can be written as 



S c (x,y) 



%(T 



c 



i\)(x)i\)(y) 



(2.13) 



From Eqs. (|2.5| ) and (|2.9| ), one can easily derive the equations of motion for the 
quark Green functions: 



(i $ x - m )S c (x, y) = S c {x, y) 



(2.14) 



i T 



c 



(o(x)T a + 7r°(x)r + 7r+(x)r_ + ir-(x)T+) ^(x)^(y) 



where the contour delta function is defined as 

{5 A (x — y) if both x and y are on the upper branch of C 
—5 4 (x — y) if both x and y are on the lower branch of C 
elsewhere 

Defining the contour ordered self energy via 

d 4 yJ: c (x,y)S c (y,z) = 

'a{x)Y a + 7r°(x)r + 7r+(x)r_ + tt~{x)T + ) ^{x)^{y) 



(2.15) 



c 



(2.16) 



-i T c 



•■c 



where the subscript C at the integral sign indicates time integration along the con- 
tour, one can cast Eq. (|2.14j) into the Schwinger-Dyson like form 



(ift x -m )S c (x, z) = 5 c (x } z) + / d yTr (x,y)S c (y, z) 

J c 

Analogously, one can derive the adjoint equation 



S c (x,z)(-z0 z -m o ) = 5 c (x } z)+ / d y S (x,y)E (y,z) 

•J c 



(2.17) 



(2.18) 



By defining retarded and advanced Green functions and self energies via 

S R (x,y) = S + (x,y) - S<(x,y) = S>(x,y) - S~(x,y) (2.19a) 

S A (x,y) = S + (x,y) - S>(x,y) = S<(x,y) - S~(x,y) (2.19b) 

£ R (x,?/) = Z + (x,y) - E<(x,y) = E>(x,y) - Z-(x,y) (2.19c) 

£ A (x,?/) = E+(x,y) - E>(x,y) = E<(x,y) - Tr{x,y) (2.19d) 

and disentangling the contour integrations in Eqs. ( |2.17| ), ( j2.1S| ), one arrives at the 
equations of motion for the Green functions 

(S, 1 ® S R > A ) (x, y) = 5\x -y)+ (e* a ® S R > A ) (x, y) (2.20a) 

(S R ' A ® S^) (x, y) = 5\x -y)+ (S R > A ® (x, (2.20b) 

(S^ 1 <g> 5 > ' < ) (x, y) = (S > ' < ® S A + S fi (g) (x, y) (2.20c) 

S > ' < ®S b 1 )(x,y) = (S > ^ < ®J: A + S R ®^ > ' < )(x,y) . (2.20d) 



In order to abbreviate the notation, the convolution operator ®, which is defined by 

(A ® 5) (x, *) = / A z) (2.21) 

and the inverse of the bare Green function 

S b \x, y) = (i & - m )6\x - y) (2.22) 

are introduced in Eqs. ( |2.20|) . Note that 'inverse' here means inverse with respect 
to the convolution fl2.21| ): 

(F&F- 1 ) (x,y) = 5\x-y) . (2.23) 

The general framework for deriving transport equations from Eqs. ( |2.20| ) has been 
outlined in Refs. Q], an explicit calculation for the NJL model has been given in 
Refs. |T3|Jl^| . Since the aim of the present work is the derivation of transport equa- 



tions for mesons, these calculations are not repeated here. It should however be 
mentioned that with the definition of the self energy (|2.16 ), £ (a;, y) still contains 



a Hartree part, which in the case of the NJL model leads to the spontaneous break- 
down of chiral symmetry and the appearance of constituent quarks. For further 
details, the reader is referred to Refs. [|T0| , |T3f| . In Ref. [14]], a gradient expansion of 



Eqs. ( |2.20| ) is performed and subsequently the quasiparticle ansatz for S 



>,< 
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S < (x,p) 



S > (x,p) 



in 5ff>5cc' 



(p' + m q (x)) 5(p - E q (x,p))n q (x,p) 



E q (x,p) 2N c N f 

-5 (p + E g {x,p)) (2N c N f - n q (x, -p)) 



(2.24a) 



in 5ff>5cd 



(jtt+m q (x)) 5 (p - E q (x,p)) {2N c N f - n q {x,p)) 

(2.24b) 



E q (x,p) 2N c N f 

-5 (p + E q (x } p)) n q (x, -p) 



with E q (x,p) = Jp 1 + m? q (x), is employed in order to derive the Boltzmann equation 



d t + d p E q (x,p)d x - d x E q (x,p)dp n q (x,p) 

dp 



(2.25) 



f .^Xr [S < (x, p)S > (x,p) — S > (x,p)S < (x,p)] 
Jo 2t[ 



for the quark density n q (x,p). The integration on the right hand side serves to 
eliminate the delta functions contained implicitly in S ,>,< . The dynamical quark 
mass m q (x) in Eq. (|2.25| ) has in turn to be calculated from the gap equation 

d 3 p 2N c Nj — n q (x, p) — n q (x, p) 



m q (x) = m + 2Gm q (x) 



(2.26) 



/|p1<a (2vr) 3 E q (x,p) 

Together with an expansion of the nondiagonal selfenergies S >,< , to be discussed 
later, Eqs. ( |2.25|) and ( |2.26|) provide a set of equations for the computation of n q (x,p) 
and n q (x,p). A numerical evaluation of these equations has been given in Refs. |[15|| . 



B. Meson Fields 



In analogy to Eq. ( |2.13| ) the Green functions for the meson fields are defined via 



A£ fry) = -i{(T c W{x)a{y)\) - (a(x)) (a(y))} 
A%(x,y) = -i {(T c [^(x)7r (y)]) - (n°(x)) (n°(y))} 
A c ± (x,y) = -i {(Tc [n+(x)n-(y)}) - (tt+(x)) (^(y))} 



(2.27a) 
(2.27b) 
(2.27c) 



for sigmas, neutral pions and charged pions, respectively. Employing the equations 

), one 



of motion for the fields derived from ( |2.6j ) and the commutation relations ( pT8 



obtains the equation of motion for the contour ordered mesonic Green functions 



2G 



ed 2 Xo + ^-)A°(x,y) = 8 c (x,y) 



(2.28a) 



i{(T c [^(x)r^(x)a(y)]) - (^(x)T^(x)) (a(y))} 
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i {(T c [$(x)r i/>{x)Tr°(y)]) - ($(x)Toi>(x)) (rc°(y))} 



- [£d 2 X0 + — )A < i(x,y)=5 c (x,y) 



(2.28b) 



f2.28c) 



-i {(T c [^(x)T + ^(x)n-(y)}) - ($(x)T + ij>(x)) (tt"^))} . 

Defining the contour ordered polarization function U^(x,y) for the sigma field via 

/ d t yHS(x,y)ASfaz)= (2.29) 
J c 

-i {(t c [^(x)T a ^(x)a(z)]) - (^(x)T^(x)) (a(z))} , 

and analogously the polarization functions U.q(x, y) for neutral pions and H±(x,y) 
for charged pions, Eqs. ( |2.28| ) can be rewitten as 



1 

2G 



^l + ^)A c k (x,z) 



S c (x, z) + [ d*y n%(x, y)A c k {y, z) (2.30) 



for k = a, 0, ±. The adjoint equation is 
1 

2G 



6 c (x,z)+ f d i yA < jj(x,y)I?j!(y,z) . (2.31) 
Jc 

In the limit e — > 0, the time derivative on the left hand side of Eqs. (|2.30|) and fl2.31|) 



drops out and one reobtains the equations of motion of Ref. [17 



Disentangling the contour integration and introducing retarded and advanced 
quantities as in Eq. ( |2 . 1 9| ) leads to 



A,- 1 ® A? A ) (x, y) = 6\x -y) + (nf A ® A? A ) (x, y) 
(A? A ® A," 1 ) y) = 5\x -y)+ (A*' A ® Of - A ) (*, y) 



Ar 1 ® A^'< 



A 



(x,y) 
fay) 



>,< 



A A k 



nf ®A k 



>,< 



>,< 



fay) 



Tlt + A^Up < )(x,y) , 



where the inverse of the bare propagator is given by 

1 



At\x,y) 



2G 



5(x - y) 



(2.32a) 
(2.32b) 
(2.32c) 
(2.32d) 



(2.33) 



Equations Q2.32| ) are the mesonic analogue of Eqs. ( |2.20| ) and will form the basis 
of the following investigations. A major difference between Eqs. ( [2.32| ) and the 
corresponding equations for particles, which occur as elementary degrees of freedom 
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in the microscopic interaction model, is that the inverse bare propagator defined in 
Eq. (|2.33|) does no longer contain derivatives, as does e. g. the inverse bare quark 
propagator given in Eq. ( |2.22|) . It will be seen in the following, however, that this 
does not inhibit the description of mesons within the framework of transport theory. 



III. TRANSPORT EQUATIONS FOR MESONS 



The derivation of transport equations from Eqs. ( |2.32|) can in principle be done 
in the standard fashion outlined in Refs. PUT"^. The drawback of this approach is, 
however, that the quasiparticle ansatz, which for bosons reads 

in 



Ek(x,p) 



ITi 



E k (x,p) 



n k (x,p)5(p - E k (x,p}) 
+ (1 +n- k (x, -p))5(p + E k (x,p)) 
(1 + n k (x,p))6(p - E k (x,p)) 
+n k (x, -p)S(p + E k (x,p}) 



(3.1a) 



(3.1b) 



does not solve the equations of motion ( p.32c|) , fl2.32d| ) unless one neglects off-shell 
terms ||. Since the problem treated here differs from the standard formalism in con- 
sidering the propagation of non-elementary particles, it is not clear a priori, whether 
this procedure is justified. It is thus desirable to use a formalism which consistently 
eliminates off-shell terms within the framework of the gradient expansion, as has 
been developed in Ref. ||. It has been shown there, that the unwanted terms can 
be eliminated by considering off-shell corrections to the ansatz (|3.1|) . Thus, in the 
following the formalism developed in || will be applied to ( [2.32| ). 



A. Mesonic Spectral Functions 

In order to transform Eqs. ( |2.32| ) to momentum space, a Wigner transformation 
is performed, which for a function F(x, y) of two coordinate arguments is defined 
by 

F(x,p) = J d 4 ue ipu F(x + u/2,x-u/2) . (3.2) 

The Wigner transform of a convolution of two functions can be evaluated employing 
a gradient expansion Q . Keeping only terms up to first order gradients, one obtains 
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d 4 ue ipu (F®G)(x + u/2,x-u/2) (3.3) 

= F(x,p)G(x,p) + -{F(x,p);G(x,p)} , 

introducing the Poisson bracket of two functions in momentum space, 

{F(x,p)\ G(x,p)} = d p F(x,p)d x G(x,p) - d x F(x,p)d p G(x,p) . (3.4) 

In the following, also the Wigner transform of triple and quadruple convolutions will 
be needed, which can be easily inferred from Eq. ( |3.3|) . 

By transforming Eqs. ( |2.32a| ) and ( p.32b|) and taking the sum and the difference 
of the transformed equations, one obtains 

- ^ A k A (x,P) = 1 + n«> A (x,p)A*> A (x,p) (3.5a) 

= {nf A (x,p);Af A (x,p)} . (3.5b) 

Equations ( |3.5|) are solved by 

A% A (x,p) = ^L- . (3.6) 

l + 2GU^ A {x,p) 



The mesonic quasiparticle energy and the width can be identified as the location of 
the poles of the retarded Green function (|3.6|) , i. e. as the solution of the dispersion 
relation 



l + 2Gnf(a;,p) = . (3.7) 



In the vicinity of a pole, the denominator of Eq. (|3.6j ) can be Taylor-expanded with 
respect to po as 

1 + 2GU%{x,p) « - 2G . = 2E k (x,p) (po - E k (x,p) + ^T k (x,pj) , 

(3.8) 

defining the quasiparticle energy E k (x,p), the quasiparticle width T k (x,p) and the 
effective meson-quark coupling g k (x,p) for meson species k. With Eq. ( |3.8|) , A k (x,p) 
becomes in the vicinity of the pole 

A«(*,P)« 9l(x,p)-ia k (x,p) _ (3Q) 

2E k (x,p) [po - E k (x,p) + |r fc (x,p) 
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and the spectral function becomes 

Pk (x,p) = i (A*(x,p) - Af(x,p)) (3.10) 

_ 1 \T k (x,p)g 2 k (x,p) + (go - E k (x,p))a k (x,p) 
E k (x,p) (po - E k {x,p)) 2 + |r^(x,p) 

Note that is a real function, since in momentum space A^ is the complex conju- 
gate of A^, as can be inferred from its definition [[|. 

It is easy to see that the imaginary part of the pole residue, a k , is of the same 
order as the width IV Thus, in the limit 1^ — > one obtains p k ~ 5(po — E k ), 
i. e. the mesons are strictly on shell. The expansions ( |3.9| ) and ( |3.10| ) contain only 
a particle pole at positive p , whereas one should expect also an antiparticle pole at 



negative po from symmetry reasons |L8[ . Since these poles are separated by a finite 
gap, however, the form given above will be sufficient. 

For illustration, the spectral function for the pion is shown in Figs. and |3] 
for thermal equilibrium and p = at T = and T = 300 MeV, respectively, 
in the random phase approximation [10]. In this case, the spectral function is a 



function of po only. One can distinguish two cases: at low temperatures as in 
Fig. 0, Eq. (|3.7|) has solutions for real \po\ < 2m q . In this case the spectral function 
has the form of a delta function in the vicinity of the solution of Eq. ( |3.7| ), i.e. 
contains contributions of bound states. Furthermore, one observes the appearance 
of a continuum at |po| > 2m g , which stems from quark-ant iquark scattering states. 

At high temperatures, Eq. (|3.7| ) has no longer solutions for real, but only for 
complex po with |3?(po)| > 2m q . Physically this means that the pion becomes 
unstable at high temperatures. Details about this so-called Mott transition can be 
found in Ref. |HJ. The spectral function at a temperature of T = 300 MeV is 
shown in Fig. |3|. One notices, that the pion does not manifest itself by a delta peak 
any longer, but rather by a resonance peak. The dashed line in Fig. |] shows the 
pole approximation according to the second line of Eq. ( |3.10j ), which can be seen 
to give a good approximation to the full spectral function even at this rather high 
temperature. Note that it is important to include these resonant states in the theory 
since they give a contribution to the effective number of degrees of freedom in the 



thermodynamical quantities, as has been shown in Ref. p4 |. 
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B. Mesonic Densities 



The determination of A^' A allows one to extract information about the properties 
of mesons. These functions do not, however, give direct information about the 
particle densities. These must be extracted from the off-diagonal Green functions 
Aj^. As was stated above, the approach developed in Ref. M will be applied. This 
formalism starts by returning to the equations of motion in coordinate space and 
observing that Eqs. (|2.32cj ) , ( |2.32d| ) are equivalent to || 



A? 



n 



>,< 



A^ 



(3.11] 



where the coordinate arguments have been dropped for simplicity. Analogously, p k 
can be rewritten as 



Pk 



*(a£-a£ 



i{A 



a,<; 



Af ® 7fc ® A^ 



(3.12) 



In Eq. ( |3.12j ), the auxiliary quantities 



n fc 



^(nf + m 1 



2 

• (n 



• (n 



n fe <; 



(3.13a) 
(3.13b) 



have been introduced. Since in momentum space the retarded and advanced quanti- 
ties are complex conjugate to each other [|]], Eq. ( |3.13| ) is for this case commensurate 
with a decomposition of H k ' A into its real and imaginary part. 

Starting from the identities (|3.11|) and (|3.12p , it is possible to decompose A 
and pk into a singular and a regular part 



>,< 
k 



Pk — Pk,s + Pk, 



(3.14a) 
(3.14b) 



which in turn are defined by 



A>,< 



AJ;< 



Pk,s 
Pk,r 



^(A^®n>'<® PA; -p fc 



n 



A^ 



k 



Af + A 



A 



Ik® Pk ~ Pk® Ik 



A£ ® lk <g> Af + A£ ® lk ® 



n 

A 



.4 



(3.15a) 
(3.15b) 
(3.15c) 
(3.15d) 
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The physical implications of the decompositions ( |3.14 ) and ( p. 15 ) have been 
discussed in Ref. ||. The decomposition of A fc ' and p k into a singular and regular 
part corresponds to a separation of long time and short time phenomena. It is thus 
possible to describe the singular part (which corresonds to the long time phenomena) 
in terms of quasiparticles. This picture is further confirmed by observing that by 
Wigner transforming Eq. ( |3.15c|) and computing the gradients using Eq. (|3.6| ), one 
obtains p ktS = \lkp\ i n momentum space, i. e. for a given width the singular part 
of the spectral function p k}S forms a much narrower peak around the quasiparticle 
pole than the spectral function itself. 

If a description of mesons in terms of quasiparticles is desired, one thus needs 
an equation of motion for the singular part of A^ ,<: only. Following Ref. ||, this is 
achieved by computing 



A< 



in two ways: Inserting ^A 



R,A 
k 



-1 



k,s 



A," 1 



A fc < 



Ilf A from Eqs. ( gggafl , ( gggg ) leads to 



A^ 



(3.16) 



(A*)' 1 ® A< s - A< s ® (A^)" 1 ] = i ( 7fc ® A< s + A 



*L ® Ik 



(3.17) 



A,; 1 - n. 



A," 1 



n. 



On the other hand, inserting Eq. ( p. 15a ) and Eq. ( |3.12| ) into ( |3.16 ) yields 



—i 



A?\ 1 ®A < 



/,.., - A^ s ® [A k 
1 



-i 



(n<® Pfc + Pfc ®n<) 



Af ®n<® A^ + Af®n<® Af ® 7fe ) . 
Equating ( |3.17| ) and ( |3.18| ) and performing a Wigner transformation gives 

{a,- 1 - n fc; A< s } - - { lk - p k IL< (Af + At ) } = i (n<A> s - n>A< 



(3.18) 



(3.19) 



where now all functions depend on x and p. In deriving Eq. ( |3.19|) , it has been used 
that 



Pk,s(x,p) = Pk(x,p) - p k , r (x,p) 

= Pk(x,p) - -lk(x,p) (Af (x,p) 2 + At(x,pf 



(3.20) 



Equation (|3.19|) already contains the ingredients of the Boltzmann equation. A term 
similar to the first one on the left hand side of Eq. (|3.19|) appears also in the standard 
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formalism, where it gives rise to a mean field term. The second term on the left 
hand side is particular to the present approach and has to be treated further. The 
right hand side has already the form of a collision term. 

In accordance with the assumption that A fc can be described in terms of quasi- 
particles, these functions are factorized as 



A kA X 'P) = - i Pk,s( X iP)fk(x,P) 

A k,s( x ,P) = -ipk,s( x ,P) I 1 + fk(x,p)] 



(3.21a) 
(3.21b) 



where f k is assumed to be a smooth function of po- It is necessary to use the same 



function f k for both Aj5 s and A^ s in order to maintain the identity % (A^ s — A 



Pfc s . The left hand side of Eq. (|3.19| ) becomes 

{a,- 1 - n fe; A< s } - i { lk - Pk u< (Ag + a£)} = -i PKs {a,- 1 - n, ; f k ) 

-if k {A,- 1 - n fc ; p k>s } - X - { lk - p k U< (Aj? + At ) } . 
The second term on the right hand side can be transformed to || 

{A,- 1 - II fc ; p M } = | Tfe ; -^ lkPk (A* + A£)| . 
With this, one can rewrite the right hand side of Eq. (|3.23| ) as 

-ip k;S {A," 1 - n fc ; f k ) - | 7fc ; X -p k (A? + Ajf) [i lk f k + III 

+ % - lkPk (a£ + Ai) { 7fc ; f k } ■ 



(3.22) 



(3.23) 



(3.24) 



The second term here is effectively of second order in the gradients and can thus be 
neglected. The third term can be rewritten to be M 



\lkPk (Af + Af) { 7fc ; f k } = p k> 



A,- 1 - U k 



{lk, fk} 



(3.25) 



so that Eq. ( |3.19| ) becomes 



iPs,k 



{A^-n k ;f k } 



a,- 1 - n A 

ik 



{ik, fk} 



Ps , k [U<(l + f k )-U>f k 



(3.26) 



Since p ktS is strongly peaked near the quasiparticle energy, it is possible to drop 
the common factor p k)S and to consider the rest of Eq. (|3.26|) on the mass shell 
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Po = E k (x,p). This corresponds to approximating the singular part of the spectral 
function by a delta function: 



Ps,k(x,p) 



[6(po - E k (x,p)) - 5(po + E k (x,p))} 



(3.27) 



E k (x,p) 

This in turn allows the introduction of the quasiparticle ansatz ( |3.1|) as an ansatz 
for the singular part of A^ ,<: only: 



ITT 



E k (x,p) 



E k (x,p) 



n k (x,p)5(p - E k (x,p)) (3.28a) 
+(1 + n k (x, -p))5(po + E k (x,p}) 

(l + n k {x,p))S{p -E k {x,p)) (3.28b) 
+n k (x,-p)5(p + E k (x,p}) . 



These expressions differ from Eqs. ( |3.1| ) by the factor g%, which stems from the pole 
residue of the retarded Green functions, cf. Eq. (|3.9|). Similar factors appear also in 
the standard formalism, where they account for wave function renormalization ||. 
Equation ( |2.24| ) contains no such factor since the quark self energy in the Hartree 
approximation is momentum independent. 

From Eqs. ( |3.28| ), one obtains f k (x,p) = n k (x,p) at positive pq. Inserting this 
into the left hand side of Eq. (|3.26|) and using the expansion ( |3.8|) yields 



{A^-U;f k } 



2E k 



.9t + 4 



9k(Po ~ E k ) - -a k T k 



n k 



(3.29) 



9k 



d t + d p E k d x — d x E k d p 



where the terms containing a k have been neglected, since they are of second order 
in the width and thus beyond the validity of the approach. 

Within the same accuracy, the second term on the left hand side of Eq. ( |3.26| ) 
vanishes and one obtains 



dt + d p E k d x — d x E k d p 



n k 



■ 9l 
2E k 



[U<(l+n k )-U>n k ] 



(3.30) 



i. e. the Boltzmann equation for mesons. 

Equation ( |3.3U| ) is one of the major results of this paper. It shows, that transport 
theory is applicable to the formation of bound states and that the mean field term, 



18 



i.e. the left hand side of Eq. (|3.30|) , has precisely the same form as it has for 
elementary particles. This means that in a situation, where no free quarks are 
present, mesons can propagate freely and can also carry kinetic energy. At a first 
glance, this might be surprising, since the left hand sides of Eqs. ( |2.32| ) do not contain 
any derivatives, which usually are associated with the kinetic energy of the particles 
HJ. This is nevertheless only true for models, where the interaction appears as a 
small perturbation and the drift term has to be present already for the noninteracting 
case. Bound states, on the other hand, are nonperturbative phenomena, for which 
the presence of interactions is essential. The interaction part of the retarded and 
advanced Green functions does thus gain a greater importance than the bare part, 
giving rise to poles which are not present in the noninteracting case. These poles, in 
turn, define the kinetic energy which later appears in the drift term of the transport 
equations. 



IV. INTERACTIONS BETWEEN QUARKS AND MESONS 

In the last two sections, the Boltzmann equations for quark and meson degrees of 
freedom have been derived and an explicit form of the mean field (Vlasov) part has 
been given. It remains to evaluate the self energies in order to express the collision 
integrals also in terms of the particle densities. In this section, the lowest order 
terms for the off diagonal polarization U.^' and self energies S >,< are examined. It 
is shown that these diagrams describe (i) the decay of mesons into quark- ant iquark 
pairs and the corresponding recombination, M <-» qq, and (ii) the elastic scattering of 
quarks and antiquarks via meson exchange, qq <-> qq and qq «-> qq. It turns out that 
the latter processes stem from the off-shell part of Aj^ and that the interference 
of the processes qq — > M and qq — > qq leads to a modified s-channel amplitude. 
The inclusion of off-shell corrections to Aj^ is thus essential for a description of 
collisions beyond the leading order. 

Throughout this section it will be assumed that the quasiparticle energies can 
be parametrized as E(x,p) = Jp 2 + m 2 (x). 
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A. Collision Term for Mesons 



The lowest order term for the polarization in a 1/N C expansion is shown in Fig. f|. 
For definiteness, the neutral pion is considered in the following, i. e. = ^75X3. The 
off-diagonal polarization IIq is given by 



- (x, y) = -Tr \iT iS (x, y) iT iS (y, x 
= -Tr [iT iS< (x, y) iT iS> (y, x)} 



(4.1) 



After Wigner transformation this becomes 



-iU^(x,p) 



(4.2) 



d A pi d 4 p 2 



(2tt) 5 (p-p 1 -p 2 )Tr\iT iS < (x } pi)iT iS > (x } -p 2 )} 



(2tt) 4 (2tt) 4 

Inserting the quasiparticle form ( 2.24] ) into (O) and multiplying out the contribu- 
tions of quarks and antiquarks leads to the appearance of four terms, which corre- 
spond to the processes q — > qn°, q — > qn°, — > qqn (i. e. a quark, an antiquark and 
a pion are emitted spontaneously from the vacuum) and qq — > tt°. This calculation 
proceeds along the same lines as given in Ref. |T3|]. Since all momenta have to be 
taken on shell, the contributions of the first three processes vanish due to kinemat- 
ical constraints. The last process is allowed if > 2m q , i. e. for sufficiently high 
temperatures. Its contribution reads, dropping the coordinate argument x, 



-iU<(p) 



d 4 pi d 4 p 2 
(2tt) 4 (2tt) 4 
Mtt) 2 



X- 



(2vr) 4 5 4 (p -p x - pa)tr [75(^1 + m q )-y s (- fa + m q )} (4.3) 
-S(p 10 - E q (p 1 ))5(p 20 - E q (p 2 ))n q (p 2 )n q (p 2 ) . 



N c Nf 2E q (p 1 ) 2E q (p 2 

Using the on-shell constraints for both quarks and pions, this can be simplified to 

d 3 pi d 3 p 2 



-n <(,)=-l/ 



gl J (27r) 3 2E q (p l ) {2ixf2E q {p 2 ) 
x\M\ 2 n q (pi)n q (p 2 ) , 



{2nY5\p- Pl -p 2 ) 



(4.4) 



where 



\M[ 



2N c N f 



(4.5) 
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is the transition amplitude for the process qq — ► n. Note that the quantities on the 
right hand side itself depend on the state of the system. By a similar calculation as 
for Eq. ( fOp , IIq can be evaluated to be 

-* n » e> - - i f vAz) m m ^ ~ » - »> < 4 - 6 > 

x |A* | 2 (2iV c iV / - n q {p x )){2N c N f - n q (p 2 )) 



with |.M| 2 given by Eq. Q Q. Inserting Q and §J) into (|3T30| ) yields the full 
Boltzmann equation: 



d t + d p E n d x - d x E n d p \ n n (p) (4.7) 

1 r d 3 p 1 d 3 p 2 
^2eJ (27rf2E q {p 1 ) (2Tx)nE q (p 2 ) 

x {n q (px)n q (p 2 )[l + n n (p)} - n*(p)[2N e N f - n,(pi)][2JV c JV/ - n q (p 2 )}} . 



f d 3 Pi d 3 p 2 4 4 — — t 



This equation has exactly the form expected for a quark-meson plasma interacting 
via the decay of mesons into quarks and the recombination of quarks into mesons. 

B. Collision Term for Quarks 

The two lowest order diagrams for the quark self energy in a 1/N C expansion 
are shown in Fig. |5|. The diagram shown in Fig. [5]a corresponds to the Hartree self 



energy and does not contribute to the collision term |L4[], so that the lowest order 
contribution to the collision term stems from the diagram shown in Fig. |5|b. The 
contribution of this diagram is given by 

d^pi d 4 p 2 
(2tt) 4 (27r) z 

xiT k iS < (x, pi)ir£iAfc (x, -p 2 ) 



- s<(*,„) = Y,f ^h^i^^p - pi - n) (4. 



The gain term of Eq. (|2.25| ) thus reads 



2 Note that the transition amplitudes for qq — > 7r° and 7r° — > qq differ due to the different 
averaging factors for the incoming particles. If the transition amplitude for the process 
7T° — > qq would be used in Eq. ( |4.6| ), the blocking factors would have to be replaced with 
l-n q (p)/(2N c N f ). 
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Tr[£<(p)S>(p)] =£/ (2^4 -P2)^A>(-p 2 - 



(4.9) 



xTr 



In this expression, the off-diagonal meson Green function appears. For the 
calculation of E < , one has to take into account both the singular and the regular 
part of this function. The two contributions of A^ s and A^ r will be computed 
separately in the following. 

1 . Collision Term from the Singular Part 



For simplicity, only the contribution of the neutral pion is computed explicitly. 
The other contributions can be evaluated along the same lines and give similar re- 
sults. To compute the contribution of the singular part, one has to insert Eqs. (|2.24 ) 
and ( p.28| ) into Eq. ((4.9|). Since for S > {p) only the contribution at positive p is rel- 
evant, multiplying out Eq. ( |4.9| ) again leads to four contributions, which can be 
identified with the processes q — > qn°, qn° — > q, — > qqn° and n° — > qq. As was the 
case for the meson collision term, only the last term survives for m n > 2m q due to 
kinematical constraints. The contribution of this process reads, after a substitution 
Pi -> -pi, 



Tr[S<(p)5>(p)] s = 
d 4 p 1 d 4 p 2 



(4.10) 



/ 



x 



(2tt) 4 (2tt) 



9l 



7T 



EM) 



n^{p 2 )8{p 20 ~ K(P2)) 



xTr|7 5 r 3 

X75^3 



dff'dcc' 

2E q {p x ) N cN f 



ix + m q ) 5 (p 10 - E 9 (pi)) (2N c N f - n ? (pi)) 



7T 



E q (p) 



S(p -E q (p)) 



2E q (p) N c Nf 

(Ppi 



y+m q )8 (po - E q {p)) {2N c N f - n q {p)) 
d 3 p 2 



{2-KY2EM) (2n)*2E v (p2) 



(2n) A 5\p + Pl -p 2 ) 



x \M\nXP2) (?N c N f - n q {p)) (2N c N f - n q (n)) 



where \M.\ 2 is again given by Eq. (|4.5|) . The total contribution of the singular part 
Aq ^ to the collision term is thus given by 
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^Tr [E<(p)S>(p) - Z>(p)S<(p)] s = (4.11) 
1 f d 3 p 1 d 3 p 2 



i27r) 4 S\p + Pl -p 2 )\M\' 



2E q (p) J (2 7 r) 3 2E q (p 1 ) (2tt) 3 2E w (p 2 ) 

x Ko(P2) [2N c N f - n q (p)} [2N c N f - n ? (ft)] - n q {p)n q {p x ) [1 + n^Pa)]} 

and describes the process 7r° <-> gg. 



2. Collision Term from the Regular Part 



In order to obtain the contribution from the regular part of A k ' to the collision 
term, one has to insert 

A>;<(x,p) = llL>><(x,p) [A*(x,p) 2 + A£(x,p) 2 } = U>'<(x,p)U (A* (x,p) 2 ) (4.12) 



into Eq. (|4.9|) . The polarization n^ ,<: (a;,p) here has to be taken from Eq. ( [4.21) . This 
leads, after performing one of the integrations by means of the delta function, to 



Tr[E<(p)S>(p)} 



(4.13) 



d A pi d A p 2 d 4 p 3 
(27r) 4 (27r) 4 (27r) 4 ' 
x Tr [*r tS< (pO zr zS* (p)) Tr [ir zS> (pa) zr iS< (p 3 )] 



(2vr) V(p - Pl +p 2 - p 3 )3* (A? \p 2 - p 3 )- 



The next part of the calculation proceeds exactly like shown in Ref . |14| . Multiplying 
out the contributions of the individual quark Green functions leads to the appearance 
of eight terms, of which five can be dropped since they correspond to processes like 
q <-> qqq and <-> qqqq. The remaining three terms can be rearranged to give 



Tr [E>(p)S<(p)\ 
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r E q (p) 
d 3 pi d 3 p 2 



S( Po - E q {p)) 
d 3 p 3 



(4.14) 



(27i) 3 2E q (p 1 ) (27r) 3 2E q (p 2 ) (2ir) 3 2E q (p 3 ) 



(2vr) 5 (p + pi-p 2 - Ps) 



x - [<P(pi - Ps) + <t>(pi ~ P2)] n q {p 2 )n q {p 3 ) [2N c N f - n q ( Pl )} [2N c N f - n q {p)} 



+ Wi - Ps) + 4>{~V2 - Ps)] n q (p 2 )n q (p 3 ) [2N c N f - n q ( Pl )) [2N c N f - n q (p)] 
where <p(p) is a shorthand notation for 



<p(p) 



(2N C N 



(4.15) 
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Equation ( 4.14 ) would have a form corresponding to a gain term due to elastic quark- 
quark and quark- ant iquark scattering, if one could identify the collision kernel (f) with 
the squared amplitude of the diagrams shown in Fig. |[ In this case <f>{p\ — pz) would 
correspond to the t-channel exchanges of Fig. |6|a and d, <p(pi — P2) to the w-channel 
exchange of Fig. []b and <p{—p2 — P3) to the s-channel exchange of Fig. |^c. Note that 
the interference terms of these diagrams belong to a higher order in 1/N C || and do 
not arise from the self-energy shown in Fig. [5]b, as has been detailled in Ref. ]14| . 

The problem with this interpretation is, that a direct calculation of the squared 
transition amplitudes corresponding to the graphs shown in Fig. ^| shows, that the 
collision kernel obtained from these is not given by Eq. (|4.15|) , but rather by || 

* (p) = (2W p4 K ( H 2 • (416) 

For the t and w-channel exchanges, this makes no difference, since in these cases 
Aq(p) is real || and thus Eqs. (|4.15 ) and (|4.16|) yield the same result. A problem 
arises, however, with the s-channel, for which Eq. ( J4.15| ) gives effectively the wrong 

exchange propagator. The results of this effect can be seen from Fig. [7[ The solid 

2 

line in this figure shows the square of the full exchange propagator, , for s- 

channel exchange as a function of y/s at thermal equilibrium and T = 300 MeV, 
whereas the dashed line shows its replacement in Eq. (|4.15| ), 3?(A^ 2 ). Since the 
latter leads to a negative collision kernel at large y/s, Eq. ( |4.15| ) cannot be regarded 
valid collision kernel for s-channel scattering. 

The physical reason for this breakdown of the theory can be seen from the 
following: the full transition amplitude (|4.16 ) contains resonance exchange, i. e. 
contributions of the process qq — > 7r — > qq. In the present approach, however, these 
contributions have to be taken out from the elastic scattering amplitude, since they 
are already counted in the collision terms for qq «-> tt. Using the full transition 
amplitude ( |4.16| ) is thus double counting. 

It is thus to be expected, that the s-channel elastic scattering amplitude has 



to be modified near a resonance. Nevertheless, Eq. ( 4.15 ) cannot be the correct 
form since it (i) gives a negative collision kernel, (ii) modifies the collision kernel 
at momenta far off from the resonance, where no modification is to be expected 
and (iii) gives also the wrong results at low temperatures, where the production of 
pions via recombination is kinematically forbidden and thus no interference can take 



24 



place. 



The form ( f4.16| ) for the collision kernel can be obtained by setting ||T8| , p5" 



A 



>,< 



n 



>,< 



Afe 



(4.17) 



This ansatz, however, inhibits the appearance of mesons as individual degrees of 
freedom. A similar result was obtained in Ref. 



18 , where it was shown that a 



quasiparticle ansatz for the meson propagator leads to collision terms containing 
processes of the type qq <-> n, but no elastic scattering of quarks, whereas the 
ansatz Q4.17 ) leads to the elastic scattering amplitude ( |4.16| ), eliminates however 
the appearance of free mesons. Note that this problem is not a generic problem of 
a model containing mesons as bound states, but rather a general problem of theo- 
ries containing interacting quarks and mesons, as can be concluded from Ref. ||18|| . 
Nevertheless, elastic scattering processes have to be obtained from the diagram [|b. 
They cannot be derived from self energy graphs like the one shown in Fig. [| which 
would lead to the appearance of elastic scattering processes when evaluated in a 
quasiparticle approximation: This graph is already implicitly included in the graph 
shown in Fig. [j|b via the meson line and its explicit appearance would thus be double 
counting. 



The conclusion from this is, that considering off-shell corrections to A k ' is 
essential in order to obtain collision terms containing elastic scattering at all, the 
form of the off-shell corrections given in Ref. || is, however, insufficient as far as 
the s-channel scattering amplitude is concerned. 

In order to find a more suitable form of the off-shell corrections, one has to 
consider a different decomposition of A^ into a singular and regular part as that 



of Eqs. ( 3.15| ). An attempt to do this, which however leads to the appearance of 
new terms in the derivation of the Boltzmann equation, is shown in the appendix. 
This decomposition leads to the alternative form 



A^ 



Af 



(4.18) 



for the regular part of A k 
A£ p 0r,p) = 



>.< 



The auxiliary quantity A k p in Eq. ( |4.18| ) is defined by 



2E k (x,p) (po - E k (x,p) + ^T k (x,p) 



(4.19) 
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This form of the off-shell corrections improves ( (4.12 ) in so far, as the effective ex- 
change propagator is modified only in a small region around the meson mass. If one 
has Aj? p m Af , one reobtains the form ( |4.15| ) for the collision kernel. Far off from 
the pole, Ajf can be neglected and one obtains the collision kernel Q4.16Q , which in 
this region is dominated by scattering states. Equation ( ^4 . 1 8| ) has thus the property 
to interpolate between these two forms. It is nevertheless not a complete solution 
to the problem, as can be seen from the dotted line in Fig. ^, which shows that 
Eq. ([4.18 ) gives a good approximation to Eq. ( 4.16Q far off the pion mass. This is 
especially important in order to obtain the correct scattering amplitude at low tem- 
peratures, where the formation of pions by recombination is kinematically forbidden. 
Near the meson mass the transition amplitude is modified, as expected. Equation 
( |4.18| ) gives, however, still a negative contribution at this energy and can thus not 
be regarded as a complete solution to the problem. It is nevertheless worthwhile to 
look for improved off-shell corrections by considering alternative decompositions. 



V. SUMMARY AND CONCLUSIONS 



In the preceding sections, it is shown that it is possible to extend the formalism 
of transport theory in order to describe the evolution of a system of quarks and 
bound state mesons. Each of these components evolves according to a Boltzmann 
equation containing a mean field term and collision integrals, which describe the 
interactions with the other component. Since mesons appear as effective degrees of 
freedom, the dynamical transition from a quark phase to a hadronical phase can 
be modeled. The derivation of the transport equations leads to medium dependent 
transition amplitudes in the collision integrals. 

In order to obtain collision integrals for higher order processes, the consideration 
of off-shell corrections to the off-diagonal Green functions is essential. These off-shell 
corrections have to be chosen carefully, since they modify the collision kernels in a 
nontrivial way. 

In order to extend the current approach to a nonequilibrium description of the 
chiral phase transition, further work has still do be done. Although the present 
approximation allows for the transformation of quark pairs into hadrons via the 
process qq — > ir, this process cannot be expected to be very efficient, since its 
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probability is limited by the phase space. One has thus to consider the processes 
qq — > 7T7T || or qq — > tittti |26| , which can proceed via a much less limited kinematics. 
This in turn demands the inclusion of higher order self energies [0]. For this work, 
it might also be necessary to introduce off-shell corrections for the quark Green 
functions. 

One of the requirements for this expansion is that it has to be symmetry con- 
serving. This concerns especially chiral symmetry and the validity of the Goldstone 
theorem. A symmetry conserving expansion for the equilibrium case has been out- 
lined in Ref . . A generalization of this work to non-equilibrium and the derivation 
of collision terms for this expansion is left for a future publication. 
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APPENDIX: DERIVATION OF IMPROVED OFF-SHELL CORRECTIONS 

The pole approximation of the retarded Green function given in Eq. ( |3.9| ) reads 

Ajfc.p)* fx,p)-ia k (x,f) (A1) 

2E k (x,p) (po - E k (x,p) + %T k (x,p) 



The imaginary part a k (x,p) of the pole residue appearing here is of the order of the 
width and thus gives no contribution to the final result in Section [D 1B[ It has also 
not been considered in Ref. ||. In the following it will thus be attempted to derive 
transport equations using the effective form 

AL(x, P) = — 9 ^Al — (A2) 



2E k (x,p) (po - E k (x,p) + ^T k (x,p) 
for the retarded Green functions. To do this, define the auxiliary quantities 

PkA x > p) = i { A k, P ( x > p) - A fc, P ( x ' p)) ( A3a 
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n fcjP (a;,p) = - (Il ktP (x,p) + Il kp (x,p) 

Jh,p( x , p) = 1 { n k,p( x l p) - u k,p( x , p)) ■ 



(A3b) 

(A3c) 
(A3d) 



It can be easily verified, that after a transformation to coordinate space & k ' p fulfill 
the equations 



(A,- 1 ® A* A ) (x, y) = 6\x -y) + (ng 



R,A 
k,p 

R.A 



(Ag 4 ® A,- 1 ) (x, y) = 5\x -y) + (A*, 
and pfc iP obeys the equation 



A 



H,A 
fc,p 

H,A 
fc,p 



foy) 



A 



k, P ® 7fc,i 



A 



/i'.p 



(A4a) 
(A4b) 

(A5) 



The decomposition of A fc ' and p k is done such that the singular parts are defined 
as 



A>,< 



Pfc, 



2 ;A« p ®n>'< 



Pk,p Pk,p 



n 



>,< 



7fc ® Pfc, P - Pk,. 



(A6a) 
(A6b) 



and the regular parts are defined as the difference of the full function and the singular 
part. 

Analogously to Eq. ( |3.17[ ), one obtains 



(a 



R 

k,p 



A L - A< 



A t 



-1 



1 



-( 7fc)P (8)A< s + A< g 



(A7) 



- (n£ ® PKp + p fciP 



^A,- 1 - Yi Kp ) ® A< s - A< s ® (A,- 1 - n fc , 
and Eq. (|3.18|) is replaced by 

-* [{ A k, P ) 1 ® A M " A M ® ( A ty 

-- ( TfcjP (g Aj p ® n< ® A£ p + Af; p ® n< ® A£ p ® TfejP ) 
After performing the Wigner transformation, this becomes 

{a,- 1 - n fc>p ; a< s } + 1 {n< (A*, + A^) p fc , p ; 7fc , P } 

Pk ti 



(A8) 



(A9) 



n fc < 



^-,-( A " ; ; • A /.-./, 



- 7fc, P A fc,, 
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By using the explicit representation (|A2|) and the definitions (|A3|) , it is possible to 
show that 

Pk, P ~ (Kp 2 + A ip 2 ) = Pk,s = \lKpp\ v ■ (A10) 

After employing the factorization Af s = —ip ktS fk, the first term on the left hand 
side of Eq. (^) becomes 

{A,- 1 - IL kjP ; A< s } = -ip k>s {A,- 1 - IT fciP ; f k ) - if k {A,- 1 - II fciP ; p kjS ] (All) 
= -ip h>8 { A^ 1 - n fciP ; f k } - if k jA^ 1 - n fcjP ; - (j k - ^ Kp ) p\ \ 

-ifk \lk, P \ \lk, P (Ag p + Ajjp) p k , P \ , 

which generalizes Eq. ( |3.23| ) to the present case. Repeating the steps which led to 
the derivation of Eq. ( |3.26| ) leads now to 

-ipk,s {^b 1 - n fc , p ; f k } - ifk | A^ 1 - n fcjP ; - (7* - 7^) pj^J (A12) 

+ipk, s — — — {7fc, P ; fk} - {ik, P \ \pk,p (a£ p + Af p ) {ifkik, P + n^)} 

1k,p 14 J 

= ^kPk,s + Hk,pfkPk,s ■ 

The fourth term on the left hand side can be dropped using the same arguments 
as before. The second term on the left hand side is new. Since it is proportial to 
the gradient of the difference 7^ — 7^, it is also regarded as belonging to higher 
orders and thus neglected. The remaining terms contain either a factor p kiS or p ks . 
Although these are not equal, they are of the same order and become strongly peaked 
near the quasiparticle energy in the small width limit. It is thus justified to drop 
these factors and one obtains 

{A,- 1 - n fc , p ; f k ) - Ab ~ Uk ' p { 7fc , p ; f k } = iU< - lKp f k , (A13) 

1k,p 

where all momenta have to be taken on shell. The second term on the left hand side 
vanishes. Furthermore it can be inferred from Eqs. (|A2|) and (|A3|) , that up to the 
order a\ one has ^ k ,p = Ik on the mass shell. After identifying f k with the particle 
density, one reobtains the Boltzmann Equation ( p. 30 ). 

Although this derivation gives the same form of the transport equation, it differs 
in the form of the regular part of A kr . By definition, one has in coordinate space 
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= - a m = A * ® n>«< ® A * ( A14 ) 

Transforming this to momentum space and neglecting the gradient terms, which are 
of the same order as the term neglected in Eq. ( |A12| ), leads to 

= ( A fe< - IpL) > ( A15 ) 

which is equal to Eq. (|4.18[) . 
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FIGURES 



FIG. 1. The Schwinger-Keldysh contour C in the complex time plane. 



FIG. 2. Pion spectral function as a function of po at thermal equilibrium and T = 0. 
The model parameters used for this calculation are thq = 5 MeV for the current quark 
mass, GA 2 = 2.105 for the coupling constant and A = 653 MeV for the NJL cutoff, which 
gives a quark mass of m q = 313 MeV from Eq. (2.26) and a pion mass of m n = 134 MeV 



from Eq. (3/7). Note the delta peaks at po = ±m n , which correspond to bound state 
pions. The continua for |po| > 2m 9 stem from qq scattering states. 



FIG. 3. Pion spectral function as a function of po at thermal equilibrium and 
T = 300 MeV. The model parameters are the same as for Fig. |2|, leading at this tem- 
perature to m q = 21 MeV and m n = 366 MeV. The solid line corresponds to the full 
spectral function, the dashed line gives a pole approximation. The delta peaks visible in 
Fig. H have now obtained a finite width and the pion has become a resonant state. 



FIG. 4. Lowest order diagram for the polarization 11^. Solid lines correspond to quarks. 

FIG. 5. Lowest order diagrams for the quark self energy X. Solid lines correspond to 
quarks, double lines to mesons. Diagram (a) gives the Hartree self energy, diagram (b) 
the lowest order contribution to the collision term. 

FIG. 6. Lowest order diagrams for elastic scattering by meson exchange. Solid lines 
correspond to quarks, double lines to mesons. The diagrams (a) and (b) give the t and u 
channel exchange for quark-quark scattering, (c) and (d) the s and t channel exchange for 
quark-antiquark scattering. 

FIG. 7. Exchange propagators for the s-channel of elastic quark-antiquark scattering 

2 

in equilibrium at T = 300 MeV as a function of y^s. Solid line: Aq , dashed line: 
3? (A^ 2 ) , dotted line: interpolating form according to Eq. ( |4.18| ). The vertical axis is 
given in arbitrary units. 
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FIG. 8. Self energy graph leading to elastic scattering processes when evaluated in a 
quasiparticle approximation. Solid lines correspond to quarks, double lines to mesons. 
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